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ABSTRACT 

In this paper, we present a Differential Evolution (DE) method and apply it to two problems of optimal power 
flow (OPF) and the economic load dispatch (ELD) with Valve-Point effects in Power Systems. In the first case, 
the standard IEEE 30-bus network is tested and its solution is compared to the ones solved by Particle Swarm 
Optimization (PSO), Genetic Algorithm (GA), Ant Colony Optimization (ACO) methods. For the second one, 
the NPSO is tested on 13-unit, 40-unit system and validated by comparing results with classical evolutionary 
programming (CEP), improved fast evolutionary programming (IFEP), improved particle swarm optimization 
(IPSO) and efficient particle swarm optimization (EPSO) methods. The numerical results are illustrated in 
many Figures and Tables. It has shown that the proposed method is better than the others in terms of total fuel 
costs, total loss and computational times. 
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I. INTRODUCTION 

Optimal Power Flow (OPF) and Economic Load Dispatch (ELD) problems are the important 
fundamental issues in power system operation. In essence, they are the optimization problems and their main 
objective is to reduce the total generation cost of units, while satisfying constraints. Previous efforts on solving 
OPF and ELD problems have employed various mathematical programming methods and optimization 
techniques. Recently, differential evolution (DE) algorithm has been proposed and introduced [1, 2]. The 
algorithm is inspired by biological and sociological motivations and can take care of optimality on rough, 
discontinuous and multi-modal surfaces. The DE has three main advantages: it can find near optimal solution 
regardless the initial parameter values, its convergence is fast and it uses few number of control parameters. In 
addition, DE is simple in coding and easy to use. It can handle integer and discrete optimization [1, 2]. 
In this paper, the DE is proposed for solving optimal power flow (OPF) problem. The proposed method has 
been tested on the standard IEEE 30-bus test systems [17]. The obtain results from the proposed method are 
compared to those ones from PSO [16], GA [17], ACO [19] methods. Besides, DE method is also proposed for 
solving ELD problem with valve point effects. This method has tested on 13-unit and 40-unit network. The 
obtained results are compared to those from Classical Evolutionary Programming (CEP) [4], Improved Fast 
Evolutionary Programming (IFEP) [4], Improved Particle Swarm Optimization (IPSO) [6] and Efficient 
Particle Swarm Optimization (EPSO) [5] methods. 

II. OPTIMAL POWER FLOW PROBLEM 

The OPF problem can be described as an optimization (minimization) process with nonlinear objective 
function and nonlinear constraints. The general OPF problem can be expressed as 
Minimize F(x) (1) 
subject to g(x) = 0 (2) 

h(x) <0 (3) 
where F(x) the objective function, g(x) represents the equality constraints, h(x) represents the inequality 
constraints and is x is the vector of the control variables, that is those which can be varied by a control center 
operator (generated active and reactive powers, generation bus voltage magnitudes, transformers taps, etc.). 
The essence of the optimal power flow problem resides in reducing the objective function and simultaneously 
satisfying the load flow equations (equality constraints) without violating the inequality constraints. 
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The fuel cost function is given by 

F (x) = £ (a. + b.P e . + c.P*. ) (4) 

t-i 

Where, ,/V G is the number of generation including the slack bus. P G is the generated active power at bus i. a,, bt 
and Cj are the unit costs curve for i' generator. 

While minimizing the cost function, it is necessary to make sure that the generation still supplies the load 
demands plus losses in transmission lines. Usually the power flow equations are used as equality constraints [5]. 

Tap 1 r P i (v,d)-(P G - p d ) 1 

' = I ''l = o (5) 

Where active and reactive power injection at bus i are defined in the following equation 

P(V,G) = Y V.V [g .. cos 6 . + B., sin 9.. ) (6) 

Q i (V,9) = Y J v i v j (G, 7 sin fl tf + cos B j £ jj ) (7) 

The inequality constraints of the OPF reflect the limits on physical devices in power systems as well as the 
limits created to ensure system security. The most usual types of inequality constraints are upper bus voltage 
limits at generations and load buses, lower bus voltage limits at load buses, reactive power limits at generation 
buses, maximum active power limits corresponding to lower limits at some generators, maximum line loading 
limits and limits on tap setting. The inequality constraints on the problem variables considered include 
Generation constraint: Generator voltages, real power outputs and reactive power outputs are restricted by their 
upper and lower bounds as follows 

P G , min < P Ci < P Gjmax for/ =1,2, ,N G (8) 

Q Ci mia < Q GI < Q Gl mm for i = 1, 2, ,N G (9) 

V Gi min < V a . < V 6t _ m „ for i = 1, 2, , N G (10) 

Shunt VAR constraint: Shunt VAR compensations are restricted by their upper and lower bounds as follows 

Q Ci , mii , < Q Ci < G Cl , ra „ for i = 1, 2, , N c (11) 

where N G is the number of shunt compensator. 

Transformer constraint: Transformer tap settings are restricted by their upper and lower bounds as follows 

T. min < T. < T. max for i = 1, 2, ,N T (12) 

where A^j- is the number of transformer tap. 

Security constraint: Voltages at load bus are restricted by their upper and lower bounds as follows 

v L . min < v Li < v i: mix for i = 1, 2, , N L (13) 

where Nl is the number of load bus. 

III. ECONOMIC LOAD DISPATCH PROBLEM 

The same as OPF, the economic load dispatch problem can be also described as the optimization 
(minimization) process with the following objective function 

c = £f ( (j>,) (14) 

where F. ( P. ) is the fuel cost function of the i' unit and P, is the power generated by the i unit. 
Subject to power balance constraints 
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Z p, = p d + p,.„ (15) 

where P D is the system load demand and is the transmission loss. One approach to estimate losses is by 
modeling them as a function of outputs of the system generator using Kron's loss formula of (16). 

^=^1 '•• " + E P O, B l0 +B 00 (16) 

' = i 7=1 ;'=i 

where By, B i0 , B 00 are known as the losses or /^-coefficients, 
and generating capacity constrains 

/>,,.,„ <;/>,<; P hma for; = 1,2, , N (17) 

where P, ■,,„,■„ and Pi >max are the minimum and maximum power outputs of the i' h unit. 

The smooth quadratic fuel cost function without valve point loadings of the generating units are given by (4), 
where the valve-point effects are ignored. The generating units with multi-valve steam turbines exhibit a greater 
variation in the fuel-cost functions. Since the valve point results in the ripples, a cost function contains higher 
order nonlinearity. Therefore, the equation (4) should be replaced by (18) for considering the valve-point 
effects. 

The sinusoidal functions are thus added to the quadratic cost functions as follows [6] 

F. ( P. ) = a. + b.P. + c.P. 2 + \e. x sin(/. x (P. min - P.)) (18) 

where flj, b\, C\ are the fuel cost coefficients of the i' unit and e ; and /j are the fuel cost coefficients of the f unit 
with valve point effects. 




power output, MW 

Fig.l: Example cost function with 6 valves [5] 

IV. DIFFERENTIAL EVOLUTION 



A. Overview of the DE 

In 1995, Storn and Price proposed a new floating point encoded evolutionary algorithm for global 
optimization and named it differential evolution (DE) algorithm owing to a special kind of differential operator, 
which they invoked to create new off-spring from parent chromosomes instead of classical crossover or 
mutation [1]. 

Similar to GAs, DE algorithm is a population based algorithm that uses crossover, mutation and selection 
operators. The main differences between the genetic algorithm and DE algorithm are the selection process and 
the mutation scheme that makes DE self adaptive. In DE, all solutions have the same chance of being selected 
as parents. DE employs a greedy selection process that is the best new solution and its parent wins the 
competition providing significant advantage of converging performance over genetic algorithms. 
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Fig. 2: DE cycle of stages 



B. Differential Evolution Algorithm 

DE algorithm is a population based algorithm using three operators; crossover, mutation and selection. 
Several optimization parameters must also be tuned. These parameters have joined together under the common 
name control parameters. In fact, there are only three real control parameters in the algorithm, which are 
differentiation (or mutation) constant F, crossover constant CR, and size of population NP. The rest of the 
parameters are dimension of problem D that scales the difficulty of the optimization task; maximum number of 
generations (or iterations) GEN, which may serve as a stopping condition; and low and high boundary 
constraints of variables that limit the feasible area [1, 2]. The proper setting of NP is largely dependent on the 
size of the problem. Storn and Price [1] remarked that for real-world engineering problems with D control 
variables, NP=20D will probably be more than adequate, NP as small as 5D is often possible, although optimal 
solutions using NP<2D should not be expected. In [24], Storn and Price set the size of population less than the 
recommended NP=10D in many of their test tasks. In [25], it is recommended using of NP>4D. In [26], NP=5D 
is a good choice for a first try, and then increase or decrease it by discretion. So, as a rough principle, several 
tries before solving the problem may be sufficient to choose the suitable number of the individuals. 

The DE algorithm works through a simple cycle of stages, presented in Fig. 2. 

These stages can be cleared as follow: 

1. Initialization 

At the very beginning of a DE run, problem independent variables are initialized in their feasible 
numerical range. Therefore, if the jth variable of the given problem has its lower and upper bound as Xj L and Xj U 
, respectively, then the jth component of the ith population members may be initialized as 



(0) = x. + rand (0,1). (x . - x ) 



(19) 



where rand(0,l) is a uniformly distributed random number between 0 and 1. 
2. Mutation 

In each generation to change each population member X;(t), a donor vector v;(t) is created. It is the method of 
creating this donor vector, which demarcates between the various DE schemes. However, in this paper, one 
such specific mutation strategy known as DE/rand/1 is discussed. 

X 1n x NP Parameter vectors from the current population 
o newly generated parameter vector 



F(X..>y 




Contour lines of the 
cost function 



Fig. 3: Mutation operator 

To create a donor vector Vi(t) for each ith member, three parameter vectors x rl , x l2 and x l3 are chosen randomly 
from the current population and not coinciding with the current x ; . Next, a scalar number F scales the difference 
of any two of the three vectors and the scaled difference is added to the third one whence the donor vector V;(t) 
is obtained. The usual choice for F is a number between 0.4 and 1.0. So, the process for the jth component of 
each vector can be expressed as, 



,(t) + F.(x r , At) 



,(*)) 



(20) 
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3. Crossover 

To increase the diversity of the population, crossover operator is carried out in which the donor vector 
exchanges its components with those of the current member X;(t). 

Two types of crossover schemes can be used with DE technique. These are exponential crossover and binomial 
crossover. Although the exponential crossover was proposed in the original work of Storn and Price [1], the 
binomial variant was much more used in recent applications [25]. In exponential type, the crossover is 
performed on the D variables in one loop as far as it is within the CR bound. The first time a randomly picked 
number between 0 and 1 goes beyond the CR value, no crossover is performed and the remaining variables are 
left intact. In binomial type, the crossover is performed on all D variables as far as a randomly picked number 
between 0 and 1 is within the CR value. So for high values of CR, the exponential and binomial crossovers 
yield similar results. Moreover, in the case of exponential crossover one has to be aware of the fact that there is 
a small range of CR values (typically [0.9, 1]) to which the DE is sensitive. This could explain the rule of 
thumb derived for the original variant of DE. On the other hand, for the same value of CR, the exponential 
variant needs a larger value for the scaling parameter F in order to avoid premature convergence [26]. 

In this paper, binomial crossover scheme is used which is performed on all D variables and can be expressed as 

fv. (0 if rand(Q,\) < CR 

«,.,(»> = ] (21) 

I x. . (t) else 

Ui j(t) represents the child that will compete with the parent Xij(t). 

4. Selection 

To keep the population size constant over subsequent generations, the selection process is carried out to 
determine which one of the child and the parent will survive in the next generation, i.e., at time t=t+l. DE 
actually involves the Survival of the fittest principle in its selection process. The selection process can be 
expressed as, 

fu.(t) if uu.(t)) < f(x.(0) 
x,(f + l)H (22) 
[X f (0 if f(X ( (0)<f(tf f (0) 

where, f () is the function to be minimized. From Eq. (10) we noticed that: 

• If U;(t) yields a better value of the fitness function, it replaces its target X;(t) in the next generation. 

• Otherwise, Xj(t) is retained in the population. 

Hence, the population either gets better in terms of the fitness function or remains constant but never 
deteriorates. 

C. Differential Evolution Implementation 

The proposed DE-based approach has been developed and implemented using the MATLAB software. 
Several runs have been done with different values of DE key parameters such as differentiation (or mutation) 
constant F, crossover constant CR, size of population NP, and maximum number of generations GEN which is 
used here as a stopping criteria to find the optimal DE key parameters. In this paper, the following values of DE 
key parameters are selected for the optimization of power losses and voltage stability enhancement: 

F = 0.2; CR = 0.6; NP = 150; GEN = 300 

and DE key parameters for the optimization of voltage deviations are selected as: 
F = 0.2; CR = 0.6; NP = 50; GEN = 300 

The first step in the algorithm is creating an initial population. All the independent variables which include 
generator voltages, transformer tap settings and shunt VAR compensations have to be generated according to 
Eq. (19), where each independent parameter of each individual in the population is assigned a value inside its 
given feasible region. This creates parent vectors of independent variables for the first generation. 

After, finding the independent variables, dependent variables will be found from a load flow solution. These 
ependent variables include generators reactive power, voltages at load buses and transmission line loadings. It 
should be mentioned that, the real power settings of the generators are taken from [4]. 
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V. NUMERICAL RESULTS 

A. Case of DE problem 

To verify the feasibility of the proposed DE method, the standard IEEE 30-bus system [17] has been used to 
test the OPF problem. The system line and bus data are given in [20]. The system has six generators at buses 1, 
2, 5, 8, 11, 13 and four transformers with off-nominal tap ratio in lines 6-9, 6-10, 4-12, and 28-27. The cost 
coefficients in (4) are given in TABLE. I.. The obtained results of the DE are compared with those of other 
methods as in TABLE. II. 



TABLE I. 

GENERATOR COST COEFFICIENTS OF IEEE 30-BUS SYSTEM 



Bus 


a 


b 


c 


1 


0 


2.00 


0.00375 


2 


0 


1.75 


0.01750 


5 


0 


1.00 


0.06250 


8 


0 


3.25 


0.00834 


11 


0 


3.00 


0.02500 


13 


0 


3.00 


0.02500 



820 



815 



o 810 



■ GA 

■ PSO 

■ DE 



0 50 100 150 200 250 300 

Iterations 

Fig. 4: Convergence nature of GA, PSO and DE 
in IEEE 30-bus system 

1.15 r 




0.95 - 



0 5 10 15 20 25 30 

Fig. 5: System voltage profile of IEEE 30-bus system 

The obtained results for the IEEE 30-bus system using the DE are given in Table II and the results are 
compared with those from PSO [16], GA [17] and ACO [19]. From the compared results in TABLE. II, It shows 
that the DE has succeeded in finding a global optimal solution. As visualized from the Fig.4 and Table III, it 
gives that the proposed DE method of optimization is more efficient when compared with other optimization 
methods. The optimum active power is in their secure values and is far from the min and max limits. It is also 
clear from the optimum solution that the DE easily prevent the violation of all the active constraints. The 
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security constraints are also checked for voltage magnitudes in Fig. 5. The voltage magnitudes are between their 
minimum and maximum values. No load bus was at the lower limit of the voltage magnitudes. 

TABLE H. 

COMPARISON OF RESULTS OF IEEE-30 BUS SYSTEM 



Optimization Methods 



Pi 


176.96 


179.367 


177.863 


177.1567 


Pi 


48.98 


44.24 


43.8366 


A O tLC\f\C 

48.6905 


Ps 


21.30 


24.61 


20.8930 


ii inn 

21.3013 


Ps 


21.19 


19.90 


23.1231 


ir\ m 1 a 

20.9714 


Pu 


11.97 


10.71 


14.0255 


1 1.9314 


Pn 


12.00 


14.09 


13.1199 


12.0078 


Qi 




-3.156 




-1 1.3708 


Qi 




42.543 




1 1 f f f\f\ 

32.6600 


Qs 




26.292 




30.9043 


Qi 




22.768 




34.1242 


Qn 




29.923 




18.0076 


Qn 




32.346 




9.1585 


Qc\0 


3.35 






4.9759 


Qcn 


2.20 






4.9773 


Q c is 


1.98 






4.8419 


e<i7 


3.15 






4.2934 


QclO 


4.54 






3.8339 


Qm 


3.81 






4.9725 


Qcn 


3.98 






3.2182 


QcIA 


5.00 






4.9978 


Qcl9 


2.51 






3.0210 


tlu 


1.0702 






1.0657 


na 


0.9557 






0.9211 


"15 


0.9724 






1.0012 


n 3 6 


0.9728 






0.9728 


V, 


1.0855 






1.100 


Pldss 




9.5177 


9.4616 


8.66 


Cost (S/h) 


800.41 


803.699 


803.123 


799.194 


Time (s) 




315 


20 


14.213 



B. Case of Economic Load Dispatch problem 
1. Case of thirte en-unit system 

Consider a thirteen generators case. The cost coefficients of these generators are given in [7]. The demanded 
load PD of this problem is 1800MW. 

TABLE III. 

COMPARISON OF RESULTS OF 13-UNIT SYSTEM CONSIDERING VALVE-POINT EFFECTS 



Method 


CPU time 

(sec.) 


Mean cost 
($/h) 


Maximum cost 
($/h) 


Minimum cost 
($/h) 


CEP 


293.41 


18190.32 


18404.04 


18048.21 


IFEP 


156.81 


18127.06 


18267.42 


17994.07 


FEP 


166.43 


18200.79 


18453.82 


18018.00 


MFEP 


315.98 


18192.00 


18416.89 


18028.09 


APPSO 




18014.61 


18291.92 


17978.89 


DPSO 




18084.99 


18310.43 


17976.31 


EP-SQP 




18106.93 




17991.03 


PSOSQP 




18029.99 




17969.93 


DE 


4.45 


17968.97 


17969.02 


17968.94 



TABLE IV. 

OPTIMAL DISPATCH AND THE CORRESPONDING COST IN 13-UNIT SYSTEM 



Unit P^, P,.„„, r Generation (MW) Cost ($) 

1 0 680 538.5587405 4993.5385438 
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2 0 

3 0 

4 60 

5 60 

6 60 

7 60 

8 60 

9 60 

10 40 

11 40 

12 55 

13 55 
Total Generation 




360 


150.4425834 


1547.3385496 


360 


224.3995664 


2152.8361465 


180 


109.8665500 


1129.4760320 


180 


109.8665500 


1129.4760320 


180 


109.8665502 


1129.4760359 


180 


109.8665516 


1129.4760597 


180 


109.8665500 


1129.4760321 


ISO 


109.8665500 


1129.4760321 


120 


77.3999125 


808.6529682 


120 


40.0000017 


474.5440299 


120 


55.0000000 


607.5910000 


120 


55.0000000 


607.5910000 


& Cost DE 


1800.0000 


17968.95667 



0 20 40 60 80 100 120 140 160 180 200 
Iterations 

Fig. 6: Convergence nature of SPSO, PC_PSO, SOH_PSO and DE in tested case of 13-unit system considering valve-point effects 

TABLE.III shows the minimum, mean, maximum cost achieved by the DE algorithm in 100 runs. Obviously, 
the minimum costs acquired by the proposed methods are all lower than that obtained by CEP [4], FEP [4], 
IFEP [4], MFEP [4], APPSO [7], DPSO [7], EP-SQP [21], PSO-SQP [21]. The maximum costs of DE are lower 
than the minimum costs of other method. The standard deviation of proposed methods is also lower than other 
method. These results show that the proposed methods are feasible and indeed capable of acquiring better 
solution. The optimal dispatches of the generators are listed in TABLE IV. Also note that all outputs of 
generator are within its permissible limits. 

2. Case of forty -unit system 

To verify the feasibility of the proposed DE method, the forty-unit system [7] were tested. The input data and 
the cost coefficients for 40-generating units are given in [7]. The total demanded load P D of this problem is 
10500 MW. The results obtained from the DE are compared with those as in TABLE. V. 

TABLE V. 

COMPARISON OF RESULTS OF 40-UNIT SYSTEM CONSIDERING VALVE-POINT EFFECTS 



Method 


CPU time 
(sec.) 


Mean cost 
($/h) 


Maximum 
cost (S/h) 


Minimum cost 
($/h) 


CEP 


1955.20 


124793.48 


126902.89 


123488.29 


IFEP 


1165.70 


123382.00 


125740.63 


122624.35 


APPSO 




123985.15 


126259.11 


122044.63 


DPSO 




123647.81 


125295.98 


122159.99 


IPSO 




121699.30 


122168.11 


121495.70 


EPSO 




NA 


NA 


124577.27 


MPSO 




NA 


NA 


122252.26 


ESO 




122524.07 


123143.07 


122122.16 


PSO-LRS 




122558.46 


123461.68 


122035.79 


NPSO 




122221.37 


122995.09 


122221.37 


NPSO-LRS 




122209.32 


122981.60 


121664.43 


SOH_PSO 




121853.57 


122446.30 


121501.14 
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DE 7.47 121467.99 121773.89 121416.42 



TABLE.V shows the best solution time, maximum cost, mean cost, and minimum cost achieved by the DE 
algorithms in 100 runs. The DE required the least solution time. The minimum cost achieved by DE was the 
best. The maximum costs of DE are lower than the minimum costs that obtained by CEP [4], IFEP [4], IPSO 
[6], APPSO [7], DPSO [7], EPSO [5], MPSO [23], ESO [22], PSO-LRS [24], NPSO [24], NPSO-LRS [24], 
SPSO [9], PC_PSO [9], SOH_PSO [9]. The standard deviation of proposed methods is also lower than other 
method. 

The generation outputs and the corresponding costs of the best solution are provided in TABLE. VI. The DE has 
provided better solutions compared with other methods. We have also observed that the solutions obtained by 
DE always satisfy the equality and inequality constraints. 

TABLE VI. 

PTEVIAL DISPATCH AND THE CORRESPONDING COST IN 40-UNH SYSTEM 



Unit 


P,.,mn 




Lreneration 
(MW) 


Cost ($) 


1 


36 


1 14 


110.7998 


925.0964 


2 


36 


1 14 


110.7998 


925.0963 


3 


60 


120 


97.3999 


1190.5485 


4 


80 


190 


179.7330 


2143.5503 


5 


47 


97 


87.7999 


706.5001 


6 


68 


140 


142.7998 


1604.7428 


7 


110 


300 


259.5996 


2612.8845 


8 


135 


300 


284.5996 


2779.8366 


o 


1 3 3 


inn 
JUL) 


1QA ?OOA 


T7QQ Tim 

z /ys.zjuz 


i n 
1U 


1 in 


inn 


1 ir\ r\r\r\r\ 


ZjUZ.UojU 


11 


94 


3/3 


yj.yyyo 


1 RQ1 10/11 
1 oy3 . jWj 


12 


94 


17S 

3/3 


yj.yy L t'+ 


l yuo. i joz 


13 


125 


3\l\l 




S 1 1 n 9Q71 


14 


125 


500 


304 S 1 


5149 6989 


15 


125 


500 


394.2793 


6436.5862 


16 


125 


500 


394.2793 


6436.5862 


17 


220 


500 


489.2793 


5296.7107 


18 


220 


500 


489.2793 


5288.7651 


19 


242 


550 


511.2793 


5540.9292 


20 


242 


550 


505.0000 


5627.7512 


21 


254 


550 


523.2793 


5071.2897 


22 


254 


550 


523.2793 


5071.2896 


23 


254 


550 


523.2793 


5057.2231 


24 


254 


550 


523.2793 


5057.2230 


25 


254 


550 


523.2793 


5275.0885 


26 


254 


550 


523.2793 


5275.0885 


27 


10 


150 


11.0000 


1164.0309 


28 


10 


150 


10.0000 


1140.5240 


29 


10 


150 


10.0000 


1140.5240 


30 


47 


97 


95.0000 


823.1977 


31 


60 


190 


197.9999 


1658.9037 


32 


60 


190 


197.9999 


1658.9037 


33 


60 


190 


197.9999 


1658.9037 


34 


90 


200 


180.0000 


1841.2934 


35 


90 


200 


164.7998 


1539.8703 


36 


90 


200 


205.6835 


2091.6657 


37 


25 


1 10 


109.9995 


1220.1637 


38 


25 


1 10 
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Fig.7. Convergence nature of SPSO, PC_PSO, SOH_PSO and DE in tested case of 40-unit system considering valve-point effects 

VI. CONCLUSION 

In this paper, the differential evolution method has been presented to solve the OPF and non-smooth 
ELD problems. From the obtained results of OPF and ELD are tested to IEEE 30-bus, 13-unit and 40-unit 
networks, we can present the some following conclusions + The differential evolution method has increased 
the convergence speed of our algorithm. It also means that the iterative numbers are decreased as in Fig.4. 

The proposed algorithms have the capability to obtain better solutions than various other methods in terms of 
total costs and computational times. Therefore, this algorithm is effective and efficient solving the OPF and 
ELD problems of large-scale power systems with valve point effects and FACTS devices.With the advantages of 
DE, we can use it for calculating some problems in power systems such as the OPF of AC/DC power systems, 
the OPF of interconnected systems, the FACTS location optimization, etc. All will be done in next works. 
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